Elastocapillary instability under partial wetting conditions: bending versus buckling 
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The elastocapillary instability of a flexible plate plunged in a liquid bath is analysed theoretically. 
We show that the plate can bend due to two separate destabilizing mechanisms, when the liquid is 
partially wetting the solid. For contact angles 8 e > n/2, the capillary forces acting tangential to the 
surface are compressing the plate and can induce a classical buckling instability. However, a second 
mechanism appears due to capillary forces normal to surface. These induce a destabilizing torque 
that tends to bend the plate for any value of the contact angle 9 e > 0. We denote these mechanisms as 
"buckling" and "bending" respectively and identify the two corresponding dimensionless parameters 
that govern the elastocapillary stability. The onset of instability is determined analytically and the 
different bifurcation scenarios are worked out for experimentally relevant conditions. 



I. INTRODUCTION 

Water-walking arthropods like water striders are able 
to float, despite their density, thanks to surface tension 
forces [l]]. Their superhydrophobic legs are submitted 
to a repulsive force along the contact line where the 
three phases (liquid, vapour and solid) meet. As the legs 
are long and flexible, they deform under these capillary 
forces Q ■ Figure Q] shows an experiment performed with 
extremely long artificial legs made of a soft solid, which is 
plunged into a liquid bath. One observes an elastocapil- 
lary instability that is triggered by increasing the contact 
angle 9 e : the initially immersed solid is pushed out of the 
liquid to the free surface, whenever 9 e is sufficiently large. 

In the limiting case of complete wetting, i.e. 9 e = 0, 
a compressive force is exerted on an elastic rod initially 
immersed in a liquid when its end pierces the liquid in- 
terface Such a rod buckles when the compressive 
force exceeds Eulers critical load. Consider the case of 
a plate of thickness H much smaller than its length L 
and its width W. It is submitted to a capillary force 
equal to the water surface tension 7^v" times the perime- 
ter ~ 2W. The critical force is equal to (ir/2) 2 BW/L 2 , 
where B is the bending stiffness, which can be expressed 
as B = EI/W, where E it the Young's modulus, I the 
moment of inertia and W the width of the plate. There- 
fore, buckling occurs if the plate is longer than a critical 
length L CI = (n/2) yW] {2^lv)- It is proportional to the 
elasto-capillary length 
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which is the length scale controlling a large class of elas- 
tocapillary problems (34l2|. 

One may wonder if the instability observed in Fig. [T] 
is of the same physical nature. Indeed, one can expect a 
buckling instability if the contact angle 9 e is larger than 
7r/2. Namely, the total downward force that the reser- 
voir exerts on the solid is proportional to jlv cos 9 e , and 



FIG. 1. (a,b) Photographs of an elastomeric wire of radius 
R = 300 /im and density p — 1.02 ■ 10 3 kg/m 3 brought into 
contact with a mixture of ethanol and water. From (a) to (b), 
the contact angle is increased continuously by decreasing the 
concentration of ethyl alcohol. The wire exhibits a buckling 
instability above a critical contact angle. (c,d) Independent 
determination of the contact angle, using drops of the same 
mixture on a plane substrate made of the same elastomer as 
in (a,b). On the left, the advancing contact angle is 60° with a 
50% volumic solution of ethyl alcohol in water. On the right, 
the advancing contact angle is 95° with a 5% volumic solution 
of ethyl alcohol in water. 



hence changes from "stretching" to "compressing" when 
the contact angle exceeds ir/2. However, there is a second 
mechanism that can lead to elastic deformations. Fig- 
ure [2] compares the capillary energy of an extremely flexi- 
ble object that either remains vertical or floats on the free 
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isv ~ 7sl - 1LV > 



1SV - 1SL ~ 1LV < 



FIG. 2. Capillary energy of a partially immersed plate, in the 
absence of gravity and elasticity, (a) Two sides of the plate 
are immersed, representing an energy of 2ysL per immersed 
length, (b) When bending towards the surface, one of the 
sides is no longer wetted and part of the liquid- vapor interface 
is covered. The associated energy per length is ^ysv + Jsl — 
7lv. The energy difference between (b) and (a) equals the 
spreading parameter S = jsv — 7SL — "/lv- For the partially 
wetting case S < 0, so that state (b) is energetically more 
favorable than state (a). 



threshold of instability for different 8 e . 

Our main finding is that the elastocapillary instabil- 
ity can occur even when 8 e < 7r/2, which is the regime 
where the capillary forces arc pulling on the plate and no 
"buckling" is to be expected. Indeed, this is due to the 
bending mechanism described in Fig. [2 due to the par- 
tial wettability of the substrate. In general, the threshold 
of instability depends on two dimcnsionless parameters 
that can be associated to bending and buckling respec- 
tively. Interestingly, the bending parameter is not only 
determined by the elastocapillary length £ec, but also 
involves the characteristic size of the meniscus. 

The paper is organized as follows. We first formu- 
late the elastocapillary problem and identify the relevant 
dimensionless quantities in Sec. HH In Sec. Mil we ana- 
lytically determine the threshold of instability by linear 
analysis and numerically compute the nonlinear bifurca- 
tion diagrams. The results are interpreted in experimen- 
tal context in Sec. lIVi where we also discuss the influence 
of contact angle hysteresis. The paper concludes with a 
discussion on the distribution of capillary forces in Sec.fVl 



surface of the liquid. The free energy difference is propor- 
tional to the spreading parameter S = jsv — Isl — Ilv, 
where 7sy, jsl and jlv ar e respectively the surface ten- 
sions of the solid/vapour, solid/liquid and liquid/vapour 
interfaces. As a consequence, bending is favourable under 
partial wetting conditions, S < 0, whatever the value of 
the contact angle 9 e . This is manifestly different from the 
buckling instability, which can only occur for 9 e > tt/2. 
Indeed, the mechanism for instability is not the verti- 
cally compressing force, but is a bending induced by the 
capillary torque exerted near the contact line Q . 

In this paper, we investigate theoretically the instabil- 
ity of an elastic plate plunged in a liquid of the same 
density - to avoid buoyancy. In the light of former stud- 
ies, we wish to address different issues. What are the 
mechanisms for instability: bending, buckling or a com- 
bination of the two? What is the relevant length L for 
the instability in the situation of Fig. [T] where the plate 
can be supposed to be infinite? What are the parameters 
controlling the instability? 

To illustrate the two mechanisms of elastocapillary in- 
stability, bending and buckling, we consider a long elastic 
plate that is hanging freely under the influence of grav- 
ity. We assume the thickness to be sufficiently small to 
allow for a thin plate elastic description. The bottom of 
the plate is brought into contact with a liquid reservoir 
that partially wets the liquid, with an equilibrium con- 
tact angle 6 e . To reveal the effect of surface wettability, 
we focus on the case where both sides of the plate are 
wetted by the same angle. This is fundamentally differ- 
ent from the situation prior to piercing of a rod through 
a meniscus 0-111, fo r which one of the contact lines is 
pinned to the edge of the solid - in that case the con- 
tact angle can attain any value. Our goal is to compute 
the shape of the plate and to analytically determine the 



II. ELASTOCAPILLARY FORMULATION 

The strategy of the calculation is to separately treat 
the portion of the plate that is outside the bath and the 
meniscus region near the contact line - see Fig. [3l For 
simplicity we assume that the plate and fluid are density 
matched, or equivalently, that the bottom of the plate 
reaches only just below the surface. In practice, we find 
that the characteristic length of the plate outside the 
reservoir is significantly larger than the capillary length 
that sets the size of the meniscus. This means that we can 
consider the forces and torques induced by the meniscus 
(Sec. Ill Bp as a boundary condition for the dry part of 
the plate (Sec. Ill A|) . The dimensionless equations and 
boundary condition are then summarized in Sec. Ill CI 



A. The plate outside the reservoir 

Given that the plate is very thin we can describe the 
shape by a line that we parametrize by its angle </>(s) 
with respect to the vertical direction (Fig. [3]). We use 
a curvilinear coordinate s that has s = at the level of 
the bath, and the curvature is k = d(j)/ds. In cartesian 
coordinates, we use the parametrisation x = x{z) 1 where 
z is orientated upward and z = corresponds to the 
level of the liquid reservoir. The relation between the 
two representations is: dx/ds = suk/j and dz/ds = cos0. 
We consider a very long plate that, due to gravity, follows 
the boundary condition (j)(oo) = 0, cf. Fig. [3J 

Since away from the meniscus there are no forces ap- 
plied to the plate, the shape can be computed from the 
elastica equations 0, 01 : 

B(f>" = F z sin (f>. (2) 
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FIG. 3. Left: Photograph of an elastomeric plate partially 
bent by capillary forces. We show that such bending to a sta- 
ble, finite angle is only possible due to contact angle hystere- 
sis. The plate thickness is H — 0.8 mm, its Young modulus is 
E — 66 kPa and its density is 1.02 ■ 10 3 kg/m 3 . As not partic- 
ularly clean tap water is used, the surface tension is around 
7 = 0.05 N/m. Right: Definitions of the vertical coordinate z 
and curvilinear coordinate s. The deformation of the plate is 
characterized by the deflection x(z) or local angle <f>(s). The 
characteristic curvature of the plate is iZ 1 . Inset: sketch of 
the resultant capillary forces in the meniscus region. At the 
contact lines there are pulling forces along the liquid-vapor 
interface of magnitude 7. In addition, the hydrostatic pres- 
sure pulls or pushes on the plate depending on the level with 
respect to bath. The scale of the meniscus is l-y. 



Here, B is the bending stiffness, which reflects the inter- 
nal elastic torque (per unit width) due to a curvature 0'. 
For a thin plate of thickness H and elastic modulus E 
one finds B = EH 3 /12. F z is the vertical component of 
the force (per unit plate width) on a cross-section of the 
plate and we will find below that F x = 0. The torque 
balance then reads: 



% = -Bcf)'. 



(3) 



This internal force moment is exerted by the upper por- 
tion of the plate on the lower portion of the plate. Note 
that for the situation in Fig. [3] the curvature is negative, 
0' < 0. 

The force F z on a cross-section consists of two contri- 
butions, due to gravity along the plate and due to surface 
tension at the meniscus boundary. At a location s along 
the plate, the vertical force due to gravity is simply the 
weight below s , i.e. pgHs per unit width of the plate. 



As we assume the plate to have the same density as the 
liquid, or equivalently that the bottom dips just below 
the surface, we only take into account the portion of the 
plate that is outside the reservoir. 

The capillary forces exerted in the meniscus region can 
be obtained by the virtual work principle, as sketched in 
Fig. |4j Let us consider the left side of the plate. Moving 
the plate vertically by dz, one changes the horizontal 
position of the contact line by tan 0(0), leading to an 
increase of the liquid- vapor interface on the left of the 
plate. However, this is compensated by an equivalent 
decrease of liquid- vapor interface on the right of the plate. 
A non-vanishing effect is that the vertical displacement 
increases the length of dry plate by dz/ cos 0(0), while the 
wetted part is decreased by the same amount. Assuming 
there is no contact angle hysteresis (A9 = 0), such that 
Isv — Isl = Jcos6 e , one finds the forces due to the left 
and right side: 



pL = 7(sin0(O) -cosfle) 

z ' cos 0(0) 

p R = 7(~ sin 0(0) - cos fl e ) 

2 " cos 0(0) 



(4) 
(5) 



Here we introduced the shorthand 7 = 7lv, which will 
be employed in the remainder of the paper. Again, these 
forces are per unit width of the plate. Similarly, by mov- 
ing the plate horizontally by dx, one reduces the water 
area by dx so that: 



itf = - 7 , and Ff = <y. 



(6) 



The total horizontal force F x 
vertical force reads: 



thus vanishes while the total 



F 7 _ = Ft 



F. 



R 



27 cos 9 e 
cos 0(0) 



(7) 



These results are easily generalized to incorporate con- 
tact angle hysteresis, i.e. allowing for different contact 
angles to the left and to the right of the plate @- For 
this, one replaces 2cos# e by cosO^ +cos0/j in equation 
(|7|), where 9l and Or denote the angles on left and right. 
Note that F z can be also be obtained by considering the 
force contributions around each of the contact lines (cf. 
Sec. |HB]). 

Combining ([2]) with the F z induced by gravity and sur- 
face tension derived above, one obtains the equation for 
the plate 



= pgHs + 27 



cos 9 e \ 
cos 0(0) ) 



sine 



(8) 



To analyze this equation it is convenient to introduce a 
length scale expressing the strength of gravity with re- 
spect to the bending stiffness: 



C<7#) 



1/3 



(9) 
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FIG. 4. The resultant vertical force on the plate can be ob- 
tained from the change in surface free energies due to a virtual 
displacement dz. See text for details. 



The analysis implicitly assumes that £ x is much larger 
than H. 

We show below that this is the typical scale over which 
the plate is curved. For the elastomeric plate used in 
the experiment presented in Fig. [3J £ x is around 7 mm. 
Scaling the curvilinear coordinate as 

S=j-, and 
the torque balance ([5]) becomes 



(10) 



S + 2S - 



cos 6 e 



sin <£>. 



cos 0(0) / 

Here we introduced a dimcnsionless number 



(ii) 



(12) 



The dimensionless number So expressed in (fT2j) manifests 
the importance of the vertical surface tension forces with 
respect to the bending stiffness. In the experiment pre- 
sented in Fig. [3j £ec is around 7.5 mm so that So is of 
order unity (So ^ 0.90). 

Interestingly, So can be interpreted as the ratio of £ x , 
the "effective" length of the plate, and the clastocapillary 
length I ec- This is consistent with the general picture of 
clastocapillarity, namely that surface tension can induce 
deformations (such as buckling) , when the plate is longer 
than i EC 0, B, 03 E Ell- Below we will see that under 
partial wetting conditions, there is another dimensionless 
parameter associated to the torques induced by normal 
forces. 

Equation (fTTj) can be solved analytically when the an- 
gle of deflection is small, i.e. when sin<I>(0) ~ ^(O) and 
cos$(0) ~ 1. This limit is relevant at large distances 
from the meniscus, where the plate tends to a straight 
line (Fig. [3]), as well as for describing the onset of the 



instability fScc. IIII|) . The equation then becomes 

$" = [S + 2S cos6i e ] $, (13) 
which can be solved as: 



$(S) = $i Ai(S + 2S o cos0 e ), 



(14) 



where Ai is the Airy function. The integration constant 
$i determines the amplitude of the deflection and has 
to be solved from the boundary condition at the menis- 
cus. The second Airy function Bi(S) diverges for large 
arguments and thus does not comply with the boundary 
condition <f>(oo) = 0. Using the large S asymptotics of 
Ai, wc find 



3>(S) ~ $i 



_ 2 o3/2 

e 3 

2V^SV4- 



(15) 



The plate thus naturally tends to a vertical line. Realiz- 
ing that S = s/£ x , we indeed find that l x sets the length 
scale over which the deflection decays along the upward 
direction. 



B. The meniscus region 

The plate outside the liquid is described by a sec- 
ond order ordinary differential equation and thus requires 
two boundary conditions. A first boundary condition is 
4>(oo) = 0, which, for example, was used while deriv- 
ing (fT4l) . The second boundary condition comes from 
the torques exerted at the meniscus region. As shown 
in d3J, the internal torque experienced by the plate is 
proportional to the curvature d<j)/ds, which balances the 
external torque T e applied by the meniscus: 



75 + Te = 0. 



(16) 



This boundary condition has to be evaluated at the posi- 
tion of the upper contact line. Namely, this point marks 
the edge of the domain for ([2]), for which no normal 
forces were taken into account along the plate. In the 
paragraphs below wc assume the meniscus on the left is 
higher than that on the right, as in Fig. [3J The left and 
right contact line positions, zl and zr, are found from 
the classical meniscus solutions [13 



z L = ±f 7 [2(1 - 
z R = ±£ 7 [2(1 - sin 



(0l - <X*l))] 1/2 , (17) 
+ 0( Sfl ))] 1/2 , (18) 



where the sign (symbol ±) depends on the value of the 
contact angle 9 e with respect to tt/2 or equivalently on 
the sign of cos# e . Note that we now allow explicitly 
for different contact angles on both sides of the plate. 
The length scale of the meniscus is given by the capillary 
length 



1/2 



(19) 
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and reflects the balance between surface tension and the 
hydrostatic pressure (gravity). In the conditions of Fig.[3j 
£y is around 2.3 mm. The analysis is simplified by using 
the hierarchy of length scales: 



# « £, « i v ~ * 



EC, 



(20) 



Since l 1 is significantly smaller than i x , it is natural to 
use a different scaling for the meniscus region. To avoid 
confusion with the preceding paragraph, where we scaled 
the curvilinear coordinate S = s/£ x , we scale only the 
cartesian coordinate in the meniscus region: 



z 



(21) 



One can assume £ 1 <C Iec, which suggests that the 
length scale over which the capillary forces are assumed 
to be influential is substantially smaller than the length 
scale over which the capillary-force induced bending can 
be significant. As a consequence, one can assume that 
the plate represents a negligible curvature in the menis- 
cus region, so that the angle can be considered constant, 
0(sl) = </>(s_r) = 4>{0). This gives a simple relation be- 
tween the coordinates z = s cos 0(0). 



1. Capillary forces 

Before addressing the torques, we first specify the var- 
ious capillary forces exerted by the liquid on the solid 
plate. The detailed spatial distribution of capillary forces 
is a difficult question in itself, as addressed e.g. in |18l . [l9j . 
As will be commented in Sec. IV A| the resultant forces 
can be represented as shown in the inset of Fig. [3] First, 
there is a force per length of magnitude 7 that pulls along 
the liquid- vapor interface. Second, there is a contribution 
due to hydrostatic pressure in the liquid, which is unbal- 
anced whenever the two contact lines are at a different 
height (i.e. when zj, ^ zr). This pressure is acting nor- 
mal to the solid surface and has to be integrated between 
the two contact lines. 

Projecting the tangential force contributions, one finds 
the resultant force along the plate 



F s = —7 (cos Or, + cos Or) . 



(22) 



taken in the positive s direction. For cosO e < 0, equilib- 
rium angle e > 7r/2, this force is compressing the plate. 
Similar to the classical buckling instability, such a com- 
pressive force has a destabilizing effect. For cos0 e > it 
is stabilizing. 

The normal forces add up to 



F n = 7 (sin Or - sm9 L ) + F p 



(23) 



where F p is the unbalanced hydrostatic pressure appear- 
ing on the left of the plate. Assuming that zt, > zr, 



or equivalcntly 6l < Or, this hydrostatic pressure is ob- 
tained by integration as 



Fn 



dsp(s) 



ds pgs cos 0(0) 



"2^(4-4) cos 0(0) 



Zr 



2 ' cos 0(0) ' 



(24) 



and 



where we in the last step we used pg = 7/^ 
s = Zl 1 l cos 0(0). These equation can be further worked 
out using (TIT)) , where we take 0(s_l) = 4>(sr) — 0(0). 
Combined with (|23[) this finally gives 



F„ 



-7 (cos 0l + cos Or) tan 0(0). 



(25) 



Let us emphasize that these resultant force compo- 
nents (|22l25p can be projected in the (x, z) directions 
in order to compare to the virtual work result discussed 
in Sec. Ill Al Indeed, the projection gives F x = while 
one recovers the correct F z upon replacing 2 cos e by 
cos Or, + cos Or in (J7J. This illustrates the importance of 
the force due to the unbalanced hydrostatic pressure F p . 
Its magnitude is of order 7 [see (|2"4"|) ] and involves the ex- 
pressions that depend on the contact angles 0l and Or. 
Most importantly, only by adding F p to ~/ (sin Or — sin Or,), 
one recovers the capillary forces obtained from the virtual 
work principle [see section III Aj . Therefore F p should be 
interpreted as a capillary force. 



2. Capillary torques 

Having established the capillary forces in the meniscus 
region, we are in a position to compute the associated 
torques. Since we are interested in the boundary condi- 
tion for the plate outside the reservoir, we compute the 
torque around the highest contact line, i.e. sl- Using 
the convention that positive torques induce a rotation in 
clockwise direction, the normal forces then give a torque 



i(sl - s R ) sin Or 
{Zl - Z R )sin0 R 



cos 0(0) 



(26) 



From the construction in the inset of Fig. [3] it is clear 
that this torque is destabilizing. Namely, if we consider 
a small perturbation where the plate is slightly bent to 
the right, the meniscus on the left rises higher than the 
meniscus on the right. As a consequence the surface ten- 
sion force on the right has a larger moment arm than 
that its counterpart on the left. The induced torque on 
the plate acts in the same direction as the initial pertur- 
bution, and hence, has a destabilizing effect. 
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Similarly to (|24[) . there is a torque induced by the hy- 
drostatic pressure. This is obtained by integrating over 
the pressure, now including a moment arm sl — s: 



T p = I dsp(s)(s L - s) 



where 



SL 



ds pgs cos 0(0) (s — sl) 



1 (Z L + 2Z R )(Z L - Z R ) 2 
6 7 7 cos 2 0(0) 



(27) 



Interestingly, this torque scales as (Zl — Zr) 2 , which 
reflects the fact that both the integrated pressure and 
the moment arm are proportional to Zl — Zr. For small 
asymmetry we can thus neglect T p with respect to the 
moment induced by the force at the contact line. 

Finally, the torque induced by the tangential forces is 
strictly zero when A8 = 0, as the forces act in the same 
directions. For small hysteresis, the resultant torque is 
of order ~ jHAd, since the arm for the tangential force 
is half the thickness of the plate. Clearly, this can be 
neglected with respect to T n , for which the arm is given 
by £ 7 . To summarize, we find the external torque 



Te Tn ~\~ Tp-, 



(28) 



which for small 0(0) and small hysteresis is dominated 
by T n - 



C. Dimensionless equations 

The results of the preceding paragraphs can be sum- 
marized as follows. We found that the plate in the region 
outside the bath is governed by the length scale t x , which 
in practice is much larger than the size of the meniscus 
£-y. To separate the regimes, we use the dimensionless 
curvilinear coordinate S = s/£ x outside the bath, for 
which the shape can be solved from (fTT|) . i.e. 



&'=[S+ 25V 



cos( 



cos$(0) 



sin$. 



(29) 



This it to be complemented by a boundary condition at 
S = sl/£ x w 0, since sl is of the order of £ 7 <C l x . 
This boundary condition is most conveniently expressed 
in terms of Z = zj£ 1 . When Zl > Zr this gives 

$'(0) = 



To - 



(Z L -Z R )sm6 R (Z L +2Z R ){Z L - Z f , 



cos$(0) 



while for Zr > Zl one has 



6cos 2 $(0) 



(30) 



(Z L - Zr) sin 6 L (2Z L + Z R )(Z L 



$'(0) 
Zr) 2 



cos$(0) 



6cos 2 $(0) 



(31) 



Tn = 



*-EC 



(32) 



In the experiment presented in Fig. [31 To is of order unity 
(To ~ 0.285). The values of Zl,r are determined by the 
contact angles from (|T7|) . and the plate inclination at the 
bottom $(0). The latter parameter follows as a result of 
the calculation and can be used to identify the buckling 
instability. 

Apart from the contact angle 6 e (and the hysteresis 
A#), the problem is governed by two dimensionless pa- 
rameters that can be interpreted as a ratio of length 
scales: 



Srj 



EC 



To 



f2 ■ 



(33) 



The first of these parameters can be interpreted as the 
ability to induce buckling for the tangential capillary 
force, provided that cos 9 < 0. Consistent with the stan- 
dard view of elastocapillarity, this effect is governed by 
the elastocapillary length £ec with respect to the "ef- 
fective" length of the plate, set by £ x . The second di- 
mensionless parameter sets the strength of the bending 
induced by the torque generated in the meniscus. As it 
involves a torque, the capillary length £ 7 intervenes as 
the moment arm. 



III. BIFURCATIONS 

It is clear that the straight plate, &(S) — 0, is a so- 
lution of (|29|) that satisfies the boundary condition (|30| . 
We now analyze the stability of these solutions in terms 
of the parameters So and To, for different values of 9 e . 
Throughout this section we assume no hysteresis, i.e. 
A9 = or equivalently 9 l — 8 ] r = 8 ' e . We first perform a 
linear analysis to identify the threshold and discuss the 
regimes where bending or buckling are dominant. Subse- 
quently, we numerically compute the bifurcation diagram 
by following the various solution branches in the nonlin- 
ear regime. 



A. Instability threshold 

The threshold of instability of the straight plate is ob- 
tained by linearizing the problem for small $, as already 
done in pTj) . yielding a solution 

$(5) = $iAi(5 + 2S"ocos6» e ). (34) 

Similarly, the boundary condition (|30|) can be expanded 



$'(0) ~ -2T 



cos t/„ sin ( 



^2(1- sin e ) 



$(0), 



(35) 



7 



T 



•So 



10 



$(0) o 




FIG. 5. The threshold of stability To versus So for differ- 
ent values of the contact angle 6 e = 7r/6, 7r/4, 7r/3, 7t/2. For 
cosSe > 0, the bending threshold To increases with So since 
the vertical capillary force has a stabilizing effect. The limit- 
ing case cos 9 e = has F z — 0, for which the bending thresh- 
old is independent of So- 



FIG. 6. Typical bifurcation diagram for bending induced in- 
stability, obtained from numerical integration of (I29I30I) . So- 
lutions branches characterized by the angle at the bottom of 
the plate, $(0), upon varying the bending strength To. The 
values of So = 10. and 8 e = 7r/4 were kept fixed. 



where the expression for the meniscus rise (fT7|) was used. 
Combining (|34|) and (|35j) one obtains the equation for 
a "neutral mode", which is a solution of the deflection 
profile for arbitrary (small) perturbation amplitude $i: 

sin0 e |cosfl e | Ai' [2S cos9 e } _ 
V2(l-sin0 e ) Ai[2SoCos# e ] ~ U ' W 

Indeed, this equation provides the threshold for the 
instability in terms of the parameters T , So and 9 e . This 
can be seen, e.g. by varying one of these parameters 
while keeping the other two constant. One finds that 
the internal moment Ti dominates the external torque 
T e (stable) or vice versa (unstable), as the parameter is 
varied across the neutral condition (|36| . 

B. Bending instability: 9 e < n/2 



proportional to sin6> e , and thus becomes more influential 
for larger contact angles. Second, the instability thresh- 
old increases with So, which represents the strength of 
the tangential capillary forces. In the regime 9 e < tt/2 
or cos# e > 0, these tangential forces are pulling on the 
plate and are indeed stabilizing; hence, one requires a 
larger value for To to induce the instability. In the limit 
of large So one can expand the Airy functions, yielding 

1 /2 

the asymptotics To ~ S . Finally, for the limiting case 
where the tangential capillary forces vanish, 9 e = tt/2, 
the critical T does not depend on So and can be com- 
puted analytically as 

, , . Ai'(0) 3 1 / 3 r(2/3) , , 

T tt/2 = r-r4 = -Vr^ « 0.3645 •• • 37 

K ' ' 2Ai 2r(l/3 v ' 



We now reveal the destabilizing effect of the torque in 
the meniscus, associated to the parameter T , which tries 
to bend the plate. This mechanism is most relevant for 
< it/2, for which it turns out the only destabilizing 
mechanism: the buckling parameter So is stabilizing in 
this range as it multiplies with cos 9 e > 0. In this context 
of bending, (|36|) indeed provides the critical T beyond 
which the flat solution becomes unstable. Namely, for 
larger T the torque in the meniscus T n becomes larger 
than the internal torque of the plate % for small pertur- 
bations, hence leading to instability. 

The result of the stability analysis is shown in Fig. 
depicting the critical T versus So for several contact an- 
gles below it/ 2. One observes the following trends. First, 
upon increasing 9 e the instability is triggered at a smaller 
T . This occurs as the destabilizing normal forces are 



We now analyze the non-linear behavior of the solu- 
tions above the threshold. We numerically solve (|29l30p 
and characterize the various solutions with $(0), the 
plate inclination at the bottom. Figure [S] shows a typical 
bifurcation diagram for the bending induced instability, 
by depicting the variation of $(0) with T . One recog- 
nises a supercritical pitchfork bifurcation, with a critical 
exponent 1. The diagram is for a given value of So and 
9 e , and shows that there is a critical T beyond which the 
trivial solution $(0) = becomes unstable and bifurcates 
into two stable branches. Above threshold, the plate in- 
clination saturates at a finite angle $(0). The critical 
T , can be easily read off from Fig. [5] by drawing a line 
parallel to To axis, passing through the corresponding So 
(here 5*0 = 10) and obtaining the T from the point of 
intersection of this line (here 9 e = it /A). 



0.8- 
0.6- 



7t/4- 
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-Jl/4- 
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FIG. 7. The threshold of stability T versus So for 6» e = 2tt/3. 
The vertical force F z now compresses the wire and leads to 
buckling, even in the absence of bending (To = 0). Similar 
to the classical buckling instability, the higher order branches 
correspond to all extrema of the Airy function. 



FIG. 8. Typical bifurcation diagram for the buckling induced 
instability, obtained from numerical integration of (|29I30|) . 
Solutions branches characterized by the angle at the bottom 
of the plate, <^>(0), upon varying the buckling strengths So- 
The values of To = 0.25 and 8 e = 2ir/3 were kept fixed. 



C. Buckling instability: 8 e > 7r/2 

We now consider the case 9 e > ir/2 for which the tan- 
gential forces are compressing the plate and can lead to 
the classical buckling instability. To isolate this buckling 
from the bending we can consider the case where To <C 1 
or 9 e w 7T. According to (|36|) . the onset of buckling is as- 
sociated to the rightmost maximum of the Airy function, 
i.e. Ai (cq) = 0, which gives 



So 



2 cos 9 P 



with 



CQ 



-1.01886 • 



(38) 



This value is indicated by the closed circle in Fig. [7j It 
can be seen in the figure that for To 7^ 0, the threshold is 
lowered due to the destabilizing nature of the torque in 
the meniscus region. 

Beyond the onset, one observes a sequence of branches 
associated to the other maximima and minima of the 
Airy function, at more negative arguments. In analogy 
to the classical buckling, these correspond to the higher 
order modes. The stability threshold is obtained from 
dSHJ) which gives, for T < 1, Ai'(c„) = 0, so that 



Sn = 



2 cos 6 P 



(39) 



The higher order branches correspond to large arguments 
of the Airy function, and can be determined accurately 
from asymptotics of Ai(s): 



?(4»+i: 



2/3 



(40) 



Similar to Fig. [SI in Fig. [HI we illustrate a typical bifur- 
cation diagram for the buckling induced (6 e > n/2) in- 
stability and characterize the various solution with $(0). 




FIG. 9. Shapes of the plate with $(0) = for different buck- 
ling strengths So- From left to right: So = 2.689, S = 4.302, 
So = 5.636, So = 6.830. The values of T = 0.25 and 
9 e — 27r/3 are as in Fig. [8] and were kept fixed. 



For the first branch, as we move along $(0) — 0, on 
crossing a critical Sq, the solution becomes unstable and 
produces two additional unstable solutions. This implies 
that there is no stable solution at a finite angle. Phys- 
ically, this suggests that on slight perturbation from its 
equilibrium position, the wire (or plate) ends up on the 
surface [as in Fig. [T], corresponding to $(0) = ir/2. On 
moving further along $(0) = 0, one encounters further 
bifurcations which correspond to higher order modes be- 
coming unstable. Similarly to the previous bifurcation di- 
agram, the critical values for So are obtained from Fig.[7J 
by drawing a line parallel to So axis, passing through 
T = 0.25. Typical solutions with $(0) =0 are depicted 
in Fig. [SI As in classical buckling, the successive branches 
are separated by half a wavelength. 
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IV. EXPERIMENTAL PERSPECTIVE 
A. Influence of thickness 

In this section we would like translate the analysis in 
terms of dimensionless numbers So and To to an exper- 
imental situation. As an illustration, we consider a case 
where we fix the material properties of the liquid and the 
elastic solid, and vary the thickness H. The instability 
is then reached below a critical thickness, for which the 
bending rigidity is sufficiently weak. Alternatively, one 
may perform an experiment as sketched in Fig. [lj where 
the thickness is fixed but the contact angle is varying. 

Retracing the steps of the analysis, one finds that So 
and To are constructed from the thickness of the plate 
H and the material parameters 7, E, (pg). The bending 
stiffness B by itself is not a material parameter as it 
depends on the thickness as B = EH 3 /12. By selecting 
the properties of the liquid and the elastic solid, one fixes 
two length scales 



(■EG 



E 



and (~, = — 



1/2 



(41) 



In the experiment presented in Fig. [3j (eg is around 
56 cm and the capillary length is £ 7 — 2.3mm so that 
(egI(~i — 250. We can express To and Sq in terms of the 
length scales defined in (l4"T|) as: 



So = 



(egH* 



(egH 3 



H' 7 / 3 (eg 



= H- 5/3 e 



EG 



-2/3 



-1/3 £ 2 



(42) 



(43) 



Interestingly, the length scales (eg an d (-y appear in dif- 
ferent combinations in the two parameters. We thus an- 
ticipate that the threshold thickness H presents different 
scaling laws, depending on whether the instability is due 
to bending (To) or due to buckling (So). 

In Fig. [10] we show the critical thickness for different 
values of the contact angle. We scaled the thickness by 
assuming To ~ 0(1), which implies a characteristic thick- 
ness H: 



V - f 9 / 7 f~ 2/7 
n — k- 1 i EG . 



(44) 



In the experiment presented in Fig. [3J H is around 
0.47 mm. The upper line in Fig. [10] represents the thresh- 
old thickness for instability. The main trend of the graph 
is that the instability is more difficult to reach, i.e. re- 
quires a smaller thickness of the plate, as the contact 
angle 9 e is decreased. 

A key result, however, is that even for very small con- 
tact angles there is still an instability due to the torque 
exerted by the normal forces. This contrasts the classi- 
cal buckling picture, since for small 6 e the capillary forces 
are not compressing, but are in fact pulling on the plate. 



2.5 

H 

H 2- 
1.5 

H 

0.5- 




n/4 



FIG. 10. Experimental perspective: the threshold of instabil- 
ity by varying the plate thickness H and the contact angle 8 e , 
for Ieg/I-i = 250. The plate thickness H is normalized by ~H, 
defined in (|44p as a threshold value for intermediate 8 e . For 
typical experimental conditions, H = 0.47 mm. 



Note that in this regime of small 6 e , Tq — > 00, and there- 
fore "H, as defined in (|4"4"]) . is no longer the correct length 
scale for the threshold thickness. Rather, in this regime 
the critical thickness is obtained from, 



21 1)" 3 



Ai' 



Ai 



^2(1- sin e ) 

5/3 r. \V7 



' 2^5/3 



Ceg 



1/7 



(45) 



For 6 e -> and H/H -> 0: 



■nj\ V3 



2 ^ 



5/3 

H 



■EG 



l/r 



1/2 



which yields: 



H 
H 



\(eg 



1/21 



2/3 



(46) 



(47) 



This asymptotic form is shown as the dashed line in 
Fig. [TO] Indeed, this regime involves a different com- 
bination of (eg and i-y than T-L. 

Above 7r/2, one enters the usual buckling regime. The 
graph also reveals the higher order buckling modes. From 
the previous paragraphs, a good approximation can be 
obtained as: 



H 
H 



■EG 



1/21 



— (4n 



1) 



-2/5 



(-2cos0 e ) 3/35 



(48) 

The onset is generated by the highest of the lines. Fi- 
nally, note that there is an optimal contact angle, slightly 
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before e = it, for which the instability is most easily 
reached. The scaling near e = tt turns out 



H 

u 



O \ 3/5 / „ x 3/35 



1.4988 



-EG 



3/35 



(49) 

To summarize, the critical thickness below which the 
plate becomes unstable increases with the contact an- 
gle e , except very close to e — n, where the thick- 
ness displays a maximum. Intriguingly, the dependence 
of the characteristic thickness on the material parame- 
ters is not universal, but depends on the contact angle. 
Three regimes can be identified, involving different com- 
binations of pg, 7 and E. At very small contact angles, 
one finds from (1471): 



(50) 



HOC (p ff )- 1 /3 7 2/3 jB -l/3 ee 2/3 j 

Close to e = 7r one has ([49]) : 

H oc ( A9 )- 2 / 5 7 3/5 £~ 1/5 , 



while at intermediate contact angles one has (|44|) 

tfoc( A9 )- 5 / 14 7 9/14 £- 2/7 - (52) 

B. Influence of contact angle hysteresis 



$(0) 




ji/4- 



( 51 ) $(0) 0- 



-n/4 



-Tl/2- 





(6) 







tc/4 



TC/2 



3ti/4 



Another important experimental feature is that one 
cannot eliminate a substantial hysteresis of the contact 
angle. This means that 9l ^ Or, with typical experimen- 
tal values A8 w 0.1 (in radians). Taking this hysteresis 
into account in the model, we find a small shift in the 
threshold of instability. However, the fact that hystere- 
sis breaks the left-right symmetry of the problem has a 
much more pronounced effect on the general structure of 
the solutions and their bifurcation diagrams. 

To illustrate this, we first show two bifurcation dia- 
grams corresponding to Fig. [10] in the case without hys- 
teresis (AO — 0). Figures [TTk andfTTb are both obtained 
by varying e , for two different values of the plate thick- 
ness. The upper plot corresponds to H/% = 1.54, for 
which a single bifurcation is observed. The lower plot 
corresponds to H/H = 0.77, for which the higher order 
buckling modes are crossed. 

The effect of contact angle hysteresis is revealed in 
Fig. [12l comparing the bifurcation diagrams for AO = 
and A8 = 0.1. Clearly, the left-right symmetry of the 
problem is broken by the hysteresis, as reflected by the 
splitting of the branches for $(0) > and for $(0) < 0. 
A practical consequence is that, below threshold, the 
stable solutions actually correspond to a nonzero $(0). 
This means that the stable states do not correspond to 
a straight plate: the plate is always slightly bent, even 
before the instability occurs. This effect is clearly visible 
in Fig. [31 which corresponds to a photograph of a stable 
plate that is indeed deformed towards a small angle. 



FIG. 11. Bifurcation diagram, $(0) versus # e for Ibg/I-i = 
250. The two panels correspond to different cross-sections of 
Fig. EH namely (a) H/H = 1.54, and (b) H/H = 0.77. 



V. DISCUSSION 
A. Spatial distribution of capillary forces 

The total force exerted by the liquid on the solid is 
7cos0 e per unit contact line, as follows from a thermo- 
dynamic argument based on virtual work (Sec. HT)) . In 
this paper, we have treated these capillary forces as if 
they were perfectly localized at the contact line - see the 
inset of Fig. [3l However, this is not necessarily a repre- 
sentation of the real distribution of capillary forces, since 
the solid can be submitted to a Laplace pressure wher- 
ever the solid surface is curved [HI, [20, . In the case of 
the plate, for example, this curvature is localized at the 
bottom edges and could induce an upward force due to 
a Laplace pressure on the solid. To restore the thermo- 
dynamic resultant force, this necessarily means that ad- 
ditional downward forces must be present at the contact 
line to counteract this effect. Here we will not discuss the 
origin and nature of this Laplace pressure and we refer 
the interested reader to fl8l.l2l| ; assuming that it exists, 
how are the results derived in this paper affected? 

We consider the two-dimensional situation depicted in 
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FIG. 12. Effect on hysteresis on the bifurcation diagram, 
<E>(0) versus 6 e , in the case H/U = 1.71 and Csa/i-y = 250. 
A9 = 0.1 (thick lines) is compared to the case without hys- 
teresis (thin line). The hysteresis allows for stable solutions 
of nonzero amplitude $(0). 



Fig. [T31 where the bottom of the solid can take an arbi- 
trary shape. The immersed part of the solid is submitted 
to a distribution of pressure j s k proportional to the cur- 
vature k and to the relevant surface tension coefficient 7 S . 
The total force exerted on the solid due to this curvature 
effect is written as a contour integral: 



F K = 



j s Kndl = 



f R dt „ r - 



(53) 



where / is the curvilinear coordinate, Z is the local tan- 
gent vector and n the local normal vector. Elementary 
geometry gives dt/dl = nri, where n is the local curva- 
ture of the surface. Hence, (|53|) shows that the Laplace 
contribution to the total capillary force depends only on 
the tangent vectors at the contact line and is indepen- 
dent of the shape of the immersed solid. As mentioned 
above, this must be compensated by an additional force 
at the contact line to restore the thermodynamic result. 
Similarly, the total moment of the Laplace force can be 
written as a contour integral: 

^ „ f R dr A t „ r _ n _R 
-f s rAnndl= / 7 S dl — [j s r At\ L 

(54) 

where we have used the property dr/dl = t. As for the 
resulting force, the moment also does not depend on the 
shape of the object and is equal to the moment of a force 
7 s i that would be localized at the contact line. In other 
words, neither the total force nor the total moment ex- 
erted on the upper part of the plate depend on the dis- 
tribution of capillary forces below the surface. 

We thus conclude that all results presented in this pa- 
per are perfectly insensitive to the true spatial distribu- 
tion of the capillary forces. Reciprocally, it also implies 
that the Laplace pressure on a solid cannot be charac- 
terised using bending or buckling experiments. Instead, 




FIG. 13. Sketch of the quantities required to compute the 
Laplace pressure exerted on a solid of arbitrary shape. See 
text for details. 



one must measure the deformations of the solid surface 

EL 



B. Conclusion 

In this paper we identified two separate mechanisms 
that can lead to elastocapillary instability of a flexible 
plate partially immersed in a liquid. The tangential com- 
ponents of the capillary forces can induce buckling when- 
ever the contact angle 9 e > tt/2. By contrast, the nor- 
mal components of the capillary forces, which are propor- 
tional to 7sin# e , have a destabilizing effect for arbitrary 
9 e > 0. The underlying physical mechanism can be in- 
ferred from the inset of Fig. [3J a small perturbation of 
the plate inclination induces a longer moment arm one 
one side of the plate, such that the resultant torque is 
destabilizing. Alternatively, one may consider a free en- 
ergy argument, showing that the immersed state of the 
wire is energetically unfavorable whenever the liquid is 
partially wetting (Fig. [2|). We found that the dimension- 
less number associated to this bending mechanism does 
not only involve the elastocapillary length £ec, but also 
the capillary length i 1 . The capillary length appears as 
it sets the typical moment arm for the torque. Changing 
the solid from a plate to a thin wire of radius R <§; £ c , as 
in Fig. [TJ the length scale for the moment arm becomes 
R. By estimating the physical parameters for the wire in 
Fig. [TJ we conclude that in this particular example the 
instability is not triggered by bending or buckling indi- 
vidually, but by a combination of the two mechanisms. 

The analysis of the present paper focussed on the tran- 
sition of an immersed plate, partially wetted on both 
sides, to a state where the plate is pushed to free surface. 
This can be considered as the inverse of the "piercing" 
problem that is relevant e.g. for water striders [2|. In the 
case of piercing, the contact line remains pinned on the 
edge of the solid before entering the liquid, such that one 
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side of the solid remains completely dry. It would be in- 
teresting to see if the equivalent of the stability threshold 
([55)1 could be derived for piercing as well. 

Acknowledgements: We gratefully acknowledge 



J. Bico and B. Roman for enlightening discussions. 
J.H.S. acknowledges support from Paris-Diderot Univer- 
sity where he stayed as an invited professor during the 
preparation of the manuscript. 



[1] J. W. M. Bush and D. L. Hu, "Walking on water: bio-lo- 
comotion at the interface", Ann. Rev. Fluid Mech. 38, 
33969 (2006). 

[2] K. J. Park and H.-Y. Kim, "Bending of floating flexible 
legs", J. Fluid Mech. 610, 381-90 (2008). 

[3] N. Chakrapani, B. Wei, A. Carrillo, P. M. Ajayan, 
and R. S. Kane, "Capillarity-driven assembly of two- 
dimensional cellular carbon nanotube foams" , Proc. Natl 
Acad. Sci. USA 101, 400912 (2004). 

[4] S. Neukirch, B. Roman, B. de Gaudemaris, and J. Bico, 
"Piercing a liquid surface with an elastic rod: buckling 
under capillary forces", J. Mech. Phys. Solids 55, 121235 
(2007). 

[5] F. Chiodi, B. Roman, and J. Bico, "Piercing an interface 

with a brush: collaborative stiffening", Europhys. Lett. 

90, 44006 (2010). 
[6] J. Bico, B. Roman, L. Moulin, and A. Boudaoud, "Elas- 

tocapillary coalescence in wet hair", Nature 432, 690 

(2004). 

[7] C. Py, P. Reverdy, L. Doppler, J. Bico, B. Roman, and 
C. N. Baroud, "Capillary origami: spontaneous wrapping 
of a droplet with an elastic sheet", Phys. Rev. Lett. 98, 
156103 (2007). 

[8] A. Boudaoud, J. Bico, and R. B, "Elastocapillary coa- 
lescence: aggregation and fragmentation with a maximal 
size", Phys. Rev. E 76, 060102 (2007). 
[9] C. Py, P. Reverdy, L. Doppler, J. Bico, B. Roman, and 
C. N. Baroud, "Capillarity induced folding of elastic 
sheets", Eur. Phys. J. Spec. Top. 166, 67 (2009). 

[10] B. Roman and J. Bico, "Elasto-capillarity: deforming an 
elastic structure with a liquid droplet", J. Phys. Cond. 
Matt. 22, 493101 (2010). 

[11] J. Hure, B. Roman, and J. Bico, "Wrapping an adhesive 



sphere with an elastic sheet", J. Phys. Chem. B 106, 
174301 (2011). 

[12] J. W. van Honschoten, J. W. Berenschot, T. Ondar- 
cuhu, R. G. P. Sanders, J. Sundaram, M. Elwenspoek, 
and N. R. Tas, "Elastocapillary fabrication of three- 
dimensional microstructures" , Appl. Phys. Lett. 97, 
014103 (2010). 

[13] S. P. Timoshenko and S. Woinowsky-Krieger, Theory 
of Plates and Shells 2nd edn. (McGraw-Hill, Singapore) 
(1959). 

[14] L. D. Landau and E. M. Lifshitz, Theory of Elasticity 3rd 
edn. (Butterworth-Heinemann, Oxford) (1986). 

[15] S. Tawfick, M. De Voider, and A. J. Hart, "Structurally 
programmed capillary folding of carbon nanotube assem- 
blies", Langmuir 27, 6389 (2011). 

[16] A. L. Hazel and M. Heil, "Surface-tension-induced buck- 
ling of liquid-lined elastic tubes: a model for pulmonary 
airway closure", Proc. Roy. Soc. A 461, 1847 (2005). 

[17] P. G. de Gennes, F. Brochard-Wyart, and D. Qur, Capil- 
larity and Wetting Phenomena: Drops, Bubbles, Pearls, 
Waves (Springer, New York) (2003). 

[18] S. Das, A. Marchand, B. Andreotti, and J. H. Snoeijer, 
"Elastic deformation due to tangential capillary forces", 
Phys. Fluid. 23, 072006 (2011). 

[19] A. Marchand, J. Weijs, B. Andreotti, and J. H. Snoeijer, 
"Why is surface tension a force parallel to the interface?" , 
Am. J. Phys. (in press) (2011). 

[20] E. R. Jerison, Y. Xu, L. A. Wilen, and E. R. Dufresne, 
"Deformation of an elastic substrate by a three-phase 
contact line", Phys. Rev. Lett. 106, 186103 (2011). 

[21] S. Mora, T. Phou, J.-M. Fromental, L. M. Pismen, and 
Y. Pomeau, "Capillarity driven instability of a soft solid" , 
Phys. Rev. Lett. 105, 214301 (2011). 



